SOLUTION OF THE SU(2) MANDELSTAM CONSTRAINTS 



N. J. Watson 

Centre de Physique Theorique, C.N.R.S. - Luminy, Case 907, 
F-13288 Marseille Cedex 9, France, 
email: watson@cptsul.univ-mrs.fr 



Abstract 



It is shown how the Mandelstam constraints for an SU (2) pure lattice gauge theory with 
37V physical degrees of freedom may be solved completely in terms of 37V Wilson and 
Polyakov loop variables and TV — 1 gauge invariant discrete ±1 variables, thus enabling a 
manifestly gauge invariant formulation of the theory. 
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1. Introduction 

The idea that a gauge theory may be formulated entirely in terms of gauge invariant 
variables is an appealing one, and has a long though not particularly successful history 
00 [01 00- ^ pure gauge theory, i.e. without external sources, the gauge invariant 
quantities are formed by taking the trace of matrices Ur given by the exponential of the 
vector potential integrated around closed paths T: 

TTUr = Tr P exp ^ A^^Tdx'' (1) 

where for SU{N) A^^, i = 1,2... N'^ — 1, are the vector potentials, the T* are the group 
generator matrices and P denotes path ordering. Such variables are generally known as 
Wilson loops in the lattice regularized versions of the theory. 

One of the fundamental features of such approaches is the infinite over-completeness 
of such loop variables. This implies the existence of dependences among the loops, usually 
known as Mandelstam constraints Q. These are complicated non- linear identities among 
loops which intersect or touch at at least one point, and stem just from the basic properties 
of SU (N) matrices. The form of the Mandelstam constraints for a given SU{N) depends 
on (so that going to the large- limit is not straightforward). For N = 2, which is the 
simplest case, they may all be derived from the fundamental identity[| for arbitrary SU{2) 
matrices Ua, Up 

Tr UaUp - Tr C/„ TrUfs + Tr U^U^^ = (2) 

A simple example of a Mandelstam constraint for SU (2) Wilson loops in the lattice versions 
of the theory, in this case coming directly from eqn. (2), is shown in fig. 1. 



+ 



Fig. 1. A simple example of a Mandelstam constraint for Wilson loops in SU{2) lattice gauge theory. 

The number of independent physical degrees of freedom of an SU{N) pure gauge 

theory is determined by dividing out the gauge freedom from the overall freedom of the 

system. In a lattice formulation this number is regulated; only such lattice theories will 

^The corresponding identity for any given SU{N) may be easily derived from the Cayley-Hamilton 
theorem. 
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be considered here. In the Lagrangian formulation the total number of independent 
degrees of freedom of the space-time lattice system is therefore given by the dimension 
of the quotient space '^hnksSU{N)/ ^sites SU{N). In the Hamiltonian formulation Q, in 
which a temporal gauge = is chosen and time remains continuous, the number of 
independent degrees of freedom of the spatial lattice system at a given time is similarly 
given by dim{(8)iinks'S'6''(A^)/ (^sitcs SU{N)}. In order to formulate either type of lattice 
theory directly in terms of this number of gauge invariant variables, it is necessary to be 
able to solve completely the corresponding SU{N) Mandelstam constraints in terms of a 
complete independent set of this number of loop variables (supplemented perhaps by some 
number of gauge invariant discrete variables) so that any loop can then be expressed in 
terms of this set. 

For SU{2), the two crucial identities required to solve the Mandelstam constraints 



were given by Loll |^] (here eqns. (16) and ([19|)). However, having derived the necessary 



identities, no attempt was then made to solve them directly. Instead it was argued p||10| 
that, assuming the non-existence of any global constraints, it is sufficient to show the 
mutual independence in any finite volume on the lattice of some specially chosen set of 
Wilson loop variables, where the total number of loops in this set equals the number of 
independent degrees of freedom of the system. For example, for a two dimensional lattice 
this set was taken to consist of the one and two plaquette rectangular loops at each site. 
Having demonstrated this independence, then given that these variables altogether give 
the correct number of degrees of freedom of the system and given the identities necessary 
in principle to solve completely the Mandelstam constraints, it was argued that these 
variables should give a complete description of the reduced configuration space of gauge 
invariant quantites of the system. Further discussion was given in ref. 

In this letter, it is shown how the Mandelstam constraints for SU{2) pure lattice gauge 
theory with dim{(8)iinks'S'f^(2)/ <8)sites SU{2)} = 3J\f variables may be solved directly and 
completely in terms of 3J\f continuous Wilson and Polyakov loop variables and TV — 1 
gauge invariant discrete variables taking values ±1. This enables a fully gauge invariant 
formulation of such theories. The need for these discrete variables (which have nothing to 
do with the Z2 centre of the SU{2) group) means that the argument of Loll is incorrect, 
or, more precisely, incomplete. Some possible practical consequences are briefly discussed. 

2. Solution of the SU(2) Mandelstam constraints 
Consider a set of m entirely arbitrary SU (2) matrices, for simplicity denoted generically 
now by a, /3, 7 . . . rather than Ua,Ui3,U-y . . .. Then there exist the following four identities, 
the derivations of which from eqn. (2) are explained in ref. [^: 

Tr(a) = Tr(a~^) (3) 
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Eqn. (^) is due to the fact that the representations of the SU (2) generators are reah 

TV(a/3) = Tiia)TT{P) - Tr(a/3-^) (4) 

(eqn. (|2|) triviahy rearranged). Using eqn. (^), the trace of any product involving matrices' 
inverses can always be written in terms of traces with no inverses. 

Tr(a/37) = -Tr(/?a7) + TY(a)T¥(/37) + Tr(/?)Tr(7a) + Tr(7)Tr(a/3) (5) 
-Tr(a)Tr(/3)Tr(7) 

Using eqn. (|5|), the trace of any product of matrices in a given order can always be 
expressed in terms of the trace in any other order and products of traces of products of 
fewer matrices. 

Tr(a/375) = ^ (Ti{a)Tr{(3-fS) + Tr(/3)TV(7fc) + Tr{j)Ti{6af3) + Tr{6)Ti{af3-f) (6) 
-Tr(a/3)Tr(7)TV(5) - TV(/37)TV((5)TV(a) 
-TV(75)TV(a)Tr(/3) - TV(fe)Tr(/3)Tr(7) 
+Tr(a/3)Tr(7(^) - Tr(a7)TV(/?(5) + Tr(a(5)Tr(/?7) 
+Tr(a)Tr(/3)TV(7)Tr(5)^ 



Using eqn. (P), the trace of any product of four or more matrices can always be expressed 
in terms of products of traces of one, two and three matrices. 

The identities eqns. together enable the trace of any product of matrices 

and their inverses from the given set to be expressed in terms of the traces of products 
of one, two and three matrices, which will be denoted generically by Tr(a), Tr(a/5) and 
Tr(a/37). 

However, the Tr(a/9) and Tr(af3'y) are themselves constrained. To express these con- 
straints, it is convenient to introduce the so-called L variables]^ 

L(a) = iTr(a) (7) 

L{a,l3) = i(-Tr(a/3)+Tr(ar')) (8) 

= -iTr(a/3) + iTr(«)T¥(/5) (9) 

L(a,/3,7) = ^(-Tr(a/37) + TY(arS)+Tr(a/37-^)-Tr(arS"') (10) 
+TV(a7/3) - Tr(Q7~V) - Tr(Q7/3"^) + Tr(a7"V~^)) 

= -|Tr(a/37) + iTr(«)T¥(/37) + iTV(/3)T¥(7a) + iTr(7)Tr(«/3) (11) 
-iTr(a)Tr(/3)Tr(7) 



The definitions of tfie L variables used fiere differ slightly from those used in ref. 
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These L variables have simple geometrical interpretations: writing the SU{2) matrices in 
the form 

a = a^I + icx • a (12) 

where a = (cJi, cr2i cts) are the three Pauli matrices and the real numbers ai, z = 0, 1, 2, 3, 
satisfy Og + a • a = 1 where a = (ai, 02, as) is a 3- vector, the L variables may be written 

L{a) = ao (13) 
L{a,P) = a -(3 (14) 
L{a,l3,'j) = ax (3 -J (15) 

In this form, the properties of the L variables under permutation and inversion of their 
matrix arguments are particularly transparent e.g. L{a,(3) = L((3,a). Indeed, the identi- 
ties eqns. (|^(@)(^) expressed in terms of the L variables become simply L{a) = L{a~"^), 
L{a,(3) = —L{a,(3~^) and L{a,l3,^) = — L(/3, a, 7) respectively. Also it can immediately 
be seen that L{a, (3) and L(a, /3, 7) vanish if any of their arguments is the unit matrix I. 
It is clear that these L variables are more convenient quantities with which to work than 
the original traces. The definitions eqns. (^) (|9|) (|TT|) may easily be inverted to give the 
Tr(a), Tr(Q!/3) and Tr(a/37) in terms of the L(a), L{a,(3) and L(a,/5, 7). 

There then exist two further sets of identities, given originally in ref. [^. Firstly: 

L(a, /3, 7)' + L{a, pfL{j, 7) + L(/3, 7)^^(0, a) + ^(7, a)''L{(3, (3) (16) 
-2L(a,/5)L(/3,7)L(7,a) -L(a,a)L(/3,/3)L(7,7) = 

Using eqn. ([T^), any L{a,l3,j), and hence any Tr(a/37), can always be expressed in terms 
of the L{a,P), or the Tr(a) and Tr(a/3), up to a sign. If the angle between two vectors 



ct, f3 is written then an overall product of the lengths of the vectors in eqn. (16) may 
be factored out leaving the identity in the form 



sin 6'q,^COS 9[axf3)-y + cos 0al3 + COS Of^j + COS 9.ya — 2 COS OajjCOS OfSjCOS O^a — 1 = (17) 

This identity expresses the fact that, given the angles O^p^ 0^^, 6^a between three vectors 
a, /3, 7, the relative orientation of the three vectors is specified up to a possible "reflection" 
of, say, 7 in the plane formed by a and (3. Thus there are in general two solutions ^(ax/3)7 
and vr — ^(qx/j)^ of the quadratic eqn. (|T7|). This is illustrated in fig. 2(a). The only 
information given by the L{a, (3, 7) which is not supplied by the L{a, (3) is this relative 
reflectional orientation of the sets of vectors a, (3, 7. This information is given simply by 
the sign of each of the L(a,/3, 7), so that the solution of the equation of constraint (p^) 
may be written 

L(a,/3,7) = s(a,/3,7)|L(a,/3,7)| (18) 
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where the modulus \L{a, is given entirely by the L{a,(3) as in eqn. ( [TBI ) and the 

discrete variables s(a,/3, 7) have value +1 if ^(q,x/3)7 < ^tt and —1 if ^(qx/3)7 > so 
specifying the relative reflectional orientations. This is illustrated in figs. 2(b) and 2(c) 
which show schematically in "plan view" the angles 6^13, 0^^^ d^a and the two possible 
relative orientations of the vectors. For ^(ax/3)7 = i-^- ct, /3i 7 coplanar and L(a, /3, 7) = 
0, or if any of the a,/3, 7 is the unit matrix so that again L(a,/3, 7) = 0, s(q,/3,7) is 
undefined. The s(a,/3, 7) have the same properties under permutation and inversion of 
their matrix arguments as the L(a,/3,7). 
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(a) (b) (c) 

Fig. 2(a). The two possible orientations of 7 relative to a, /3 given the three angles Qafi,Q(i~i,Q-ia- 
Fig. 2(b). The three angles Qa^^Bp-^^d^a represented schematically in "plan view" for s(q,/3,7) = +1. 
Fig. 2(c). As 2(b) but for s(a,/3,7) = -1. 

Secondly (and lastly): 

L(a, /3)2l(7, bf + L(a, ^fL{fi, bf + L(a, bfhi^H, 7)' (19) 
-L(a, /3)2l(7, 7)L(<5, 5) - L(a, 7)2l(/?, /3)L(5, 5) - L(a, 5)2l(/3, /3)L(7, 7) 
-L(/3, ifL(a, a)L{b, b) - L{(3, bfL{a, a)L{-f, 7) - L(7, bfL{a, a)L{(3, (5) 
+2L(a, a)L{l3, 7)L(7, 5)L(5, /?) + 2L(/?, /3)L(a, 7)L(7, b)L{b, a) 
+2L(7, 7)L(a, /3)L(/3, b)L{b, a) + 2L{b, b)L{a, f3)L{p, j)Li^, a) 
-2L{a, p)L{(3, 7)L(7, b)L{b, a) - 2L{a, j)L{(3, j)L{P, b)L{b, a) 
-2L{a, f3)L{l3, b)L{j, b)L{^, a) + L(q, a)L{p, l3)L{j, j)L{b, b) = 

Again, an overall product of the lengths of the vectors in eqn. (|l^) may be factored out 
leaving the identity entirely in terms of the angles between the vectors. This identity 
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expresses the fact that, given four vectors a, (3, 7, S, one of the six angles between them, 
say 9^s, may always be expressed in terms of the other five, 6^13, Oa-y, (^a&i dpyi dps in 
this case, up to a sign ambiguity: although there are altogether four possible refiectional 
orientations for the vectors 7 and 8 in the plane formed by a and /3, only two of these 
give in general distinct values for ^^5, the other two being just overall refiections of these 
two orientations, so that there are in general two solutions of this quadratic equation of 
constraint. This is illustrated in fig. 3 which shows schematically the angles in the four 
possible relative orientations of the vectors. 




s(a,/3,7) = +1 
s(a,/3,5) = +1 



s(a,/3,7) = +1 



s(a,/3,7) = -1 



s(a,/3,7) = -1 
s{a,l3,5) = -1 



(a) 



(b) 



(c) 



(d) 



Fig. 3. The four possible orientations of 7 and 5 relative to a, (3 represented schematically in "plan view" 
given the five angles 6^13, da-,, O^s, dfi-y, dfis, represented by the unbroken lines. The angle d^s is represen- 
ted by the dotted line. 



(20) 



Writing eqn. (|19|) quadratic equation for cos O^s 

a cos^0^5 + b cos 9ys + c = 



where the coefficients a,b,c are functions of the other five angles, the fact that cos 9 is 
monotonic for < 9 < ir enables the solution to be written 



cos 6*75 = -— + s(a, /3, 7) s(a, /5, 6) 
la 



(}? - 4ac)2 



2a 



(21) 



i.e. the larger angle occurs when 7 and 6 are on opposite sides of the plane formed by 
a and (3, in which case s{a, (3, 7) s{a, (3, 6) = —1. 

We now consider the number of independent continuous variables L(a) and L(a, (3) 
and discrete variables s{a, (3, 7). With m matrices a, f3,^ . . . there are m associated vectors 
a, (3, 7. . .. The lengths of these vectors are given by the m independent variables L{a) 
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via \a\ = +(1 — L(q)^)2. Given the lengths of the vectors, the information given by the 
variables L(a, f3) is the set of angles 6^13- If some particular pair of vectors, assumed non- 
parallel and of non-zero length, is chosen to form a plane, say those labelled 1 and 2, then 
there are m — 2 angles 9ia and m — 2 angles 62a required to specify the orientation of the 
remaining ni — 2 vectors relative to 1 and 2, up to the reflectional ambiguities. Including 
the angle 9i2, this means that there are 2m — 3 independent Oa/s, and hence, given the m 
L(a), 2m — 3 independent -L(a, Thus 

number of independent L{a), L{a,f3) = 3m — 3 (22) 

Similarly, there are m — 2 independent 2, a) required to specify the reflectional orien- 
tation of the m — 2 vectors relative to the plane formed by 1 and 2. Thus 

number of independent s{a,(3,'y) = m — 2 (23) 

It should be noted that, while there are 3m continuous variables associated with the 
m SU{2) matrices, there are only 3m — 3 independent continuous variables L{a), L{a, P). 
This is due to the fact that the L{a) and L{a, (3) are invariant under any rotation of all m 
vectors a, /3, 7. . .; equivalently, the Tr(a) and Tr(a/?) are invariant under a transformation 
a — > gag~^ , g G SU{2), of all m matrices a,/?, 7.... Thus the variables L{a), L{a, f3) 
are unable to describe the three "Euler angle" variables required to specify the overall 
orientation of the vectors relative to some set of axes. However, together with the m — 2 
independent s(q,/3, 7), they are fully sufficient to describe any trace formed from the m 
matrices, since such traces are themselves invariant under such an overall transformation. 

To summarize, using the formulae in this section, any trace formed from the set of 
m matrices Q,/3,7 . . . can be fully expressed in terms of some particular set of 3m — 3 
independent continuous variables L{a), L{a,(3) (or alternatively Tr(a), Tr(a/3)) and m — 2 
independent discrete variables s{a,(3,'y). Thus the SU{2) Mandelstam constaints may be 
solved completely in terms of these variables. 

3. SU(2) lattice gauge theory 
We now consider an SU{2) pure gauge system on a d-dimensional lattice with n'^ sites 
and with periodic boundary conditions^ In the Lagrangian formulation, d is the num- 
ber of space-time dimensions, while in the Hamiltonian formulation it is the number 

^ Alternatively, there are |m(m — 1) distinct L(q,/3). Setting a — 1 and /3 = 2, there are ^(m — 2)(m — 3) 
independent equations of constraint for the sets of matrices {1,2,7,5} (the left side of eqn. ([l9| ) 
vanishes trivially if any two of the a, 13, 7, 5 are equal and is invariant under any permutation of the 
a,/3,7,5). All other such eqns. ( pj| ) not involving the matrices labelled 1,2 are then dependent on these. 
Thus, given the m L[a), there are ^m{m — 1) — |(m — 2)(m — 3) = 2m — 3 independent L{a,[3). 

*A set of arguments very similar to the following can be given for the case of free boundary conditions. 
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of space dimensions, time remaining continuous. The number 3AA of independent de- 
grees of freedom of either lattice system is given by the dimension of the quotient space 
Olinks5f/(2)/0sites SU{2) i.e. 

3J\f = dim SU{2) x (number of hnks — number of sites) (24) 
= 3{d-l)n'^ (25) 

A gauge in which the gauge redundancy is reduced to the greatest possible extent 
is given by making local gauge transformations such that as many links as possible are 
fixed to the unit matrix /. Such a configuration of fixed links defines a "maximal tree" , 
characterized by the fact that, while the maximum number of links have been fixed, it is 
impossible to have any of them form a loop. An example of a maximal tree for d = 2, 
n = 4 is shown in fig. 4. The number of unfixed links on a maximal tree for any d is 
J\f +1. The "extra" three degrees of freedom coming from the unfixed "A/'+ I'th" link are 
due to the fact that it is still possible to make a global gauge transformation, rotating the 
overall system of SU{2) matrices' associated vectors relative to some frame of reference. 
Specifying this relative orientation requires three Euler angles, so that if done this would 
account for these three degrees of freedom. 



Fig. 4. A possible maximal gauge fixing tree for a 4 x 4 lattice with periodic boundary conditions. Fixed 
links are indicated by the heavy lines. 

The gauge invariant and hence "physical" quantities are the traces of loops of link 
matrices. If (and it is emphasized that this is not what will in fact be done) the TV + 1 
unfixed links on the maximally gauge fixed lattice were represented by the m matrices 
a, /5, 7 . . ., so that m = + 1, then from eqns. (|2^ ) and (p^) 

number of independent Tr(a), Tr(a/3) = 3AA (26) 
number of independent s(a, /3, 7) = J\f — 1 (27) 

Thus it is seen that when dealing only with traces and expressing them in terms of inde- 
pendent Tr(a), Tr(a/3) and s(a,/5, 7) as described in section 2, the number of continuous 
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degrees of freedom is exactly correct. 

This demonstration that for physical quantities the number of independent Tr(a) and 
Tr{al3) is correct involves some particular maximal gauge fixing - the Tr(a), Tr(a/?) and 
s(a,/3,7) used above are certainly not gauge invariant quantities. In order to arrive at 
a fully gauge invariant formulation of the SU{2) theory, rather than taking the m = 
J\f +1 matrices a, /?, 7 . . . to be the remaining links after some maximal gauge fixing, it is 
necessary to construct from the lattice links a different set of A/" + 1 independent matrices 
without reference to any gauge fixing. This set must satisfy two criteria: firstly, the trace 
of each of the matrices must be gauge invariant i.e. the matrices must each represent a 
loop; secondly, it must be possible to construct from these matrices any other loop matrix 
by taking some appropriate product of matrices and their inverses from the set. Together, 
these two criteria imply that the matrices must each consist of products of links which start 
and end at the same (arbitrary) point on the lattice. Furthermore, the second criterion 
implies that, in d dimensions, d of these matrices must consist of products of links which 
"wrap around" the periodic lattice, one in each of the directions. When traced, these 
particular matrices give Polyakov loops. An example of a set of such matrices for d = 2, 
n = 3 (hence A/" + 1 = 10) is shown in fig. 5. 



Ui 



Uo 



Us 



U5 





















U7 



Us 



Uo 



u 
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Fig. 5. A possible set of 10 basis matrices, indicated by the heavy lines, for a 3 x 3 lattice with periodic 
boundary conditions. 



It can be seen in this example that, indeed, any other loop matrix can be constructed 
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from the matrices shown. For d = 2 it is easy to see how this particular choice of matrices 
(which is only one among an infinite number) can be extended for lattices of other sizes. 
For d = 3, 4, analogous sets can be constructed, though not so easily. 

As the matrices ah start and end at the same point on the periodic lattice, it is clear 
that the trace of any product of these matrices and their inverses is gauge invariant. In 
particular, this obviously includes the Tr(a) and Tr(a/3). Crucially, this also means that 
the discrete variables s(a,/3, 7) too are manifestly gauge invariant. For example, for the 
choice of matrices shown in fig. 5, the variable s(4, 7, 10) is given from eqn. (|l8|) by 



.(4 7 10) = ^^^'^'^"^ (28) 
""^^'^''^^ |i^(4,7,10)| 



with, using the original definition eqn. (11 



L(4,7, 10) 



+- 



+- 



+- 



(29) 



(the position of the above loops on the lattice of fig. 5 is not indicated but should be 
clear) . 

Given some SU{2) lattice gauge system, it is therefore necessary to choose some basis 
set of Af+1 independent matrices a, /5, 7 . . . satisfying the above two criteria and then from 
these to choose some set of 2J\f — 1 independent Tr(a/3) and A/" — 1 independent s{a,(3,^), 
the Af + 1 Tr(a) being immediately settled. In the example of fig. 5, two of the Tr(a) are 
the Polyakov loops while the remaining eight are single plaquette Wilson loops located 
at eight of the nine elementary plaquettes of the lattice. It is immediately clear that this 
choice of basis matrices does not give a set of variables which is invariant under translation, 
90° rotation and reflection. Because, for the discrete variables, J\f — 1 = {d — l)n'^ — 1 is 
not proportional to the total number n'^ of lattice sites, it is in fact impossible to construct 
a set of variables which has the same spatial symmetries as the lattice. 

4. Conclusions 

It has been shown here how an SU (2) lattice gauge theory may be fully described in terms 
of the correct number 3M of gauge invariant continuous loop variables together with M—l 
discrete ±1 variables. Apart from the question of quantization, two practical possibilities 
immediately arise. The first concerns Monte Carlo simulations in the Lagrangian formu- 
lation. Instead of writing the action in terms of the individual link variables, it could be 
written in terms of a set of gauge invariant variables. The algebraic expression for S would 
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be extremely complicated, but the fact that it would (naturally) involve discrete ±1 vari- 
ables might enable efficient updating, similar to the extremely efficient Wolff embedding 



algorithms used for non-linear 0{N) cj-models |12]|T^]|14]. The second possibility con- 
cerns (semi-) analytic approaches in the Hamiltonian formulation. In the Kogut-Susskind 
Hamiltonian, in Schrodinger representation the electric field operators are differential op- 
erators, so that with the wavefunction ip a function of a set of gauge invariant variables, 
the Schrodinger equation Hip = Eip is a partial differential eigenvalue equation in these 
variables. On a finite lattice, such an equation might be amenable to numerical solution; 
alternatively ■0 could be expanded as a Taylor series in these variables. This latter, but 
on an effectively infinite lattice, is very similar to the approach advocated in ref. |15|. 
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